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Abstract 



The existence of a natural ultraviolet cutoff at the Planck scale is widely 
expected. In a previous Letter, it has been proposed to model this cutoff as 
I an information density bound by utilizing suitably generalized methods from 

the mathematical theory of communication. Here, we prove the mathematical 
conjectures that were made in this Letter. 



1 Introduction 



Gedanken experiments indicate that the basic spatial notions of length, area and vol- 
^ . ume lose operational meaning at the Planck scale of 10~^^m (assuming 3+1 dimen- 

sions). The argument is that any attempt at resolving smaller scales would imply a 
^ ; momentum uncertainty large enough to cause a curvature uncertainty which would 

;h ' spoil the attempted spatial resolution. It has been widely suggested, therefore, that 

^ the uncertainty relations should contain corrections that imply a finite lower bound to 

the achievable uncertainty in localizing the position or time of events, see pQ. If there 
exists such a smallest length, area or volume in nature this should imply some form of 
an upper bound to the information density that can be imprinted on physical fields. 
This is indicated as well by studies of holography, see e.g. |2|. 

The question arises, however, in which way such an upper bound to information 
density could be modelled mathematically. In [3], it was proposed to model the infor- 
mation density bound in essentially the same way that has been successful in standard 
information theory and its applications to communication engineering. 

In the early days of communication engineering, the problem arose how to quantify 
the amount of information contained in continuous signals of music, speech or images. 
Claude Shannon solved the problem in two steps in his seminal work [H] which started 
modern information theory. First, he noticed that bandlimited continuous signals such 
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as music are fully captured by discrete samples, i.e. that they are fully reconstructible 
from those samples, if the spacing of those samples does not exceed l/2a;max, where 
^max is the bandwidth of the signal, i.e., the largest frequency that occurs in the 
signal. The mathematical discipline which studies the reconstruction of functions from 
samples is now known as Shannon sampling theory and this is what we here attempt 
to generalize for applications to Planck scale physics. The proposal is that physical 
fields could be captured everywhere if sampled only at discrete points in space, as 
long as the spacing of those sampling points is sufficiently small, say of the order of 
the Planck distance. As the second step. Shannon considered the information density 
limiting effect of noise in the signal. While important, we will here not consider this 
second issue, other than mentioning that the role of noise could be played by quantum 
fluctuations. 

In order to generalize Shannon sampling theory for use in the physics of general 
curved spacetimes, it is important to develop a sampling theory that allows both the 
bandwidth and the density of samples to vary continuously: the original Shannon 
sampling theorem requires equidistantly-spaced samples but on generic curved spaces 
there are no equidistantly spaced lattices. 

In |3j, as a flrst step, a generalization of Shannon sampling theory for the case of 
one dimension was outlined and its main features were conjectured. On the basis of 
those conjectures, a new method of data compression that utilizes variable sampling 
rates (adapted to the signal) was developed and a patent has been granted, [4. Here, 
we will explicitly carry out the program which was outlined in jH] , thereby proving the 
conjectures that were made in [Hj. In addition, we show that we are dealing with a 
reproducing kernel Hilbert space and we explicitly calculate its reproducing kernel. 

Recently, in [S], this framework has been generalized to describe a covariant informat- 
ion-theoretic cutoff in arbitrary-dimensional curved spacetime. The results and tech- 
niques developed in the present paper should be very useful for explicitly carrying out 
the much more general program proposed in |B,. 

2 Terminology and notation 

We are concerned with the theory of continuous functions that can be reconstructed 
from any discrete set of samples that is sufficiently tightly spaced. This field of research, 
sampling theory, is well-established both in mathematics and in communication engi- 
neering. While we are here targeting applications to Planck scale physics we will use 
the established terminology. The translation of terms is obvious: signals mean real 
scalar fields, bandwidth means ultraviolet cutoff, the time coordinate for signals could 
be both a time coordinate or a space coordinate for fields, etc. Note that in the case 
of fields, unlike in the case of signals, samples of field amplitudes are generally not 
directly measurable. Both in quantum mechanics and in quantum field theory the 
measurement problem is deep and unsolved. If discrete samples suffice to capture a 
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field everywhere this may have profound implications on the measurement problem. 
While very interesting, we will here not further pursue this question. Also, we will be 
concerned with the sampling theory of scalar fields only. The generalization to vector, 
spinor and operator-valued fields will be considered elsewhere. 

Let us begin by noting that the mathematical language of the generalized sampling 
theory that we are concerned with is also the mathematical language of quantum 
mechanics, namely functional analysis. While this suggests the use Dirac's bra-ket 
notation, there are good reasons for using the proper mathematical notation instead. 
For example, instead of Dirac's we will write {(pjAip). The reason is that 

we will have to be precise about the domains and adjoints of operators. While an 
expression such as {A(f), ip) is unambiguous, the same expression using Dirac's notation 
would require bracketing which becomes exceedingly cumbersome in more complicated 
situations. 



3 Shannon sampling theory 

We begin by recalling elements of Shannon sampling theory. The basic Shannon sam- 
pling theorem states that in order to capture a signal (f){t) with bandwidth Umax for all 
times t, it is sufficient to record only the signals' values at the discrete points in time 
{tn} with spacing s = t^+i — tn, where: 

(1) 



If the samples are taken at this rate, which is the so-called Nyquist rate, then the 
Shannon sampling formula allows the reconstruction of the continuous signal at 
all times t from its values 0(tn) at the discrete times {tn}'- 

, ... V? sin (27r (t - Q Umax) , . 

n=—oo \ ™/ ^max 

There is a 1-parameter family of "sampling time lattices" {tn{a)} which all have the 
required spacing s, namely 

Tl 

tn{a) = ha (3) 

where a G [0,l/2u!max]- The sampling formula Eq|21 applies equally to all of these 
sampling lattices. 

The sampling theorem has a simple proof which combines discrete and continuous 
Fourier transforms. Indeed, the basic idea of the sampling theorem is reported |H| to 
go back as far as to Borel (1897), or even to Cauchy (1841). The theorem was intro- 
duced to information theory by Shannon in 1949, see jH]. In the meanwhile, numerous 
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generalizations of the theorem have been worked out. For example, in practical commu- 
nication engineering applications one usually applies a generalization of the sampling 
theorem which improves the convergence of the sampling expansion at the price of a 
slight oversampling, typically on the order of 10%. For standard references see e.g. 
[71 IH]. For recent work see, e.g., [3 El E| and references therein. The ability to 
completely recover continuous functions from discrete samples, through the sampling 
theorem, is being applied ubiquitously in fields ranging from music on CDs and sci- 
entific data taking to pure mathematics where it is used, for example, to turn certain 
sums into easier-to-evaluate integrals. 

Our aim here, as outlined in [3], is to suitably generalize Shannon sampling theory 
for application in the description of Planck scale physics in curved spacetime. To this 
end the goal is to find a sampling theorem that can accommodate variable densities of 
degrees of freedom, i.e. varying "bandwidths" and "Nyquist rates" . In the literature, 
several methods for non-equally spaced sampling are known, such as Gaussian and La- 
grangian interpolation, interlaced sampling, or Kramer's, Papoulis' or Paley- Wiener's 
sampling theorems, see, e.g., jHJin]. These, however, are either specialized or do not 
provide a direct handle on the time-varying bandwidth. Indeed, while it is clear, intu- 
itively, that the bandwidth of a class of signals can be varying in time, it is surprisingly 
difficult to give a satisfactory definition of the concept of time- varying bandwidth: 

For example, let us assume a signal obeys a constant bandlimit LVmax- If 
we consider this signal on some finite interval [ti,tf] and calculate its Fourier series 
on this interval then one might expect to find only frequency components smaller or 
equal to Umax- This is, however, not the case. Indeed, for any finite bandwidth one 
can find signals which oscillate in a given interval arbitrarily faster than the highest 
Fourier component. This shows that any attempt must fail that tries to define a 
signal's time- varying bandwidth as the maximal frequency that it contains in a moving 
temporal window. There are surprising implications for quantum fields which in this 
sense "superoscillate" , see (TU]. 

There is a further reason for why the usual notion of bandwidth does not lend itself 
to generalization to obtain a notion of time-variable bandwidth: the usual definition 
of bandwidth as the range of the Fourier spectrum allows to derive Shannon's sam- 
pling theorem, but all Fourier series - they are needed in the proof - are necessarily 
equidistant. 

For this reason, our approach here, as suggested in jH], is to aim instead for a 
generalized sampling theory based on a well-defined notion of variable Nyquist rate 
- the minimum rate at which samples must be taken can be variable. To this end, 
we will replace the use of Fourier techniques by the use of more powerful functional 
analytical techniques. Sets of sampling points will be spectra of self-adjoint operators: 
the spectra of self-adjoint operators do not have to be equidistant. 

Technically, our starting point is the observation that if is bandlimited then 
the signal t(f){t) obeys the same bandlimit. This means that no matter what bandwidth 
we choose, the corresponding class of bandlimited signals is invariant under the action 
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of the time operator T : (f){t) — > t(f){t). (This is because T is the derivative operator 
on frequency space and if a function has support only strictly within a given interval, 
then so does its derivative.) The sampling points will be the spectra of the self-adjoint 
extensions of the time operator T. 

Technically, the classes of signals with time- varying bandwidth are the domains of sim- 
ple symmetric operators T with deficiency indices (1, 1) whose self-adjoint extensions 
possess purely discrete spectra. The sampling kernel (in the special case of Eq|21 this 
is the term of the form sin(t)/t) is recognized as representing the matrix elements of 
the unitary operators which interpolate the eigenbases of the family of self-adjoint ex- 
tensions of T. The sampling kernel can be calculated explicitly in terms of the chosen 
time- varying Nyquist rate or in terms of the time- varying "bandwidth" , defined as the 
inverse of the time- varying Nyquist rate. 



4 Sampling theory in Hilbert space 

4.1 Sampling kernels as unitary operators which interpolate 
self-adjoint extensions 

The basic idea in |3J is to reduce sampling problems to the study of a certain class 
of operators, namely the class of simple symmetric operators, T, with deficiency in- 
dices (1,1) whose self-adjoint extensions T{a) have purely discrete spectra {t„(a)}. 
Here, a G [0, 27r] is a real parameter which parametrizes the ?7(l)-family of self-adjoint 
extensions T(a). We denote their eigenvalues and eigenvectors by t„(Q;) and t„(a), i.e. 

T{a)tn{a) = tn{a)tn{a). (4) 

As we will see, each such operators possess the property that the spectra {tn{a)} of 
its self-adjoint extensions are covering the real line exactly once. As a consequence, 
after fixing the phases of the eigenvectors, each vector (p in the domain of T is uniquely 
represented as a function 

cP{t) = {t,^). (5) 
This immediately provides the connection to sampling theory: 

Any vector (p is fully determined by its coefficients in one Hilbert basis, and its co- 
efficients in any arbitrary Hilbert basis can then be calculated. Here, if a Hilbert 
space vector's coefficients in the Hilbert basis {t„(a)} of one self-adjoint extension of 
T are known, then the coefficient of that vector with any arbitrary eigenvector t 
of any arbitrary self-adjoint extension of T can be calculated. Therefore, each func- 
tion 0(t) = (t, 0) is already fully determined for all t if only its values on one of the 
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spectra are given. For any arbitrary fixed choice of a, we obtain tlie reconstruction- 
from-samples formula 

0(t) = ^G'(t,t„(a))0(t„(a)), (6) 

wliere tlie sampling kernel t„(a)) reads 

G'(t,t„(«)) = (t,t„(«)). (7) 

We have merely inserted a resolution of the Hilbert space identity. The sampling ker- 
nel is therefore identified as consisting of the matrix elements of the unitaries which 
connect the eigenbases of the various self-adjoint extensions. 

In order to make this observation into a practical sampling theorem it will be neces- 
sary to establish the relation between the sampling lattices, i.e. the family of spectra 
{tn{(y)}, and the time- varying bandwidth ujmaxit)- This will be done in SecEl using 
the inverse relation between the bandwidth and the Nyquist sampling rate. Before, 
however, let us further clarify the underlying functional analytic structure. 

4.2 Simple symmetric operators 

We saw that the new and very general sampling theorem follows if we can show that 
for each simple symmetric operator T with deficiency indices (1,1) whose self-adjoint 
extensions have purely discrete spectra, the family of spectra {tn{(y)} of its self-adjoint 
extensions T{a) covers all of IR exactly once. 

To this end, let T be a closed symmetric operator on a dense domain in a 
complex Hilbert space H: 

{(j),T(f)) e M V0 G Dt, (8) 

Let us choose the scalar product to be antilinear in the first factor. Recall that T is 
called simple symmetric, see e.g. ^H], if it is not self-adjoint, and if it does not have 
any self-adjoint restriction to an invariant subspace. Note that all symmetric operators 
are closable. 

As we indicated above, the case of most interest here is that of simple symmetric 
operators T which have deficiency indices (1, 1), i.e. for which: 

dim((T ± il) Dt) = 1 (9) 

It is well known that each such operator has a ?7(l)-family of self-adjoint extensions 
T(a), where a labels the self-adjoint extensions, e.g. exp(ia) G U{1). We consider 
here only such T with self-adjoint extensions whose spectra are discrete. 

Note that for finite deficiency indices all self-adjoint extensions have the same continu- 
ous spectrum. Thus, if one of the self-adjoint extension has a purely discrete spectrum, 
so do all self-adjoint extensions. As a consequence, the spectral kernel of each such 
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operator T is empty. This is because on the one hand any continuous part in its spec- 
tral kernel would also exist in the spectrum of its self-adjoint extensions. On the other 
hand, a simple symmetric operator also cannot have any discrete part in its spectral 
kernel, since this would contradict its simplicity (T would then have self-adjoint restric- 
tions to eigenspaces) . Therefore, for T, each real number A G IR is a point of regular 
type, in the terminology of jTHj, i.e.: 

3k{\)>0 W (PeDt: ||(T- A1) 011 > fc(A) II0II (10) 

It is known that, generally, if A is a real point of regular type for a simple symmetric 
operator T with deficiency indices (n, n), then there exists a self-adjoint extension T{a) 
for which A is an eigenvalue of multiplicity n. Thus, A is also an eigenvalue of T*, of 
multiplicity at least n. On the other hand, if T is any symmetric operator with finite 
and equal deficiency indices (n, n), and if A G IR is not in its point spectrum, then the 
multiplicity of the eigenvectors of T* to the eigenvalue A does not exceed the deficiency 
index n. Thus, if A is a real point of regular type of the symmetric operator T then 
the degeneracy of the eigenvalues of T* is indeed exactly the deficiency index n. 

Here, we are considering operators with deficiency indices (1, 1) for which the entire 
real line consists of points of regular type. Thus, for each t G IR there exists a self- 
adjoint extension, T(a), for which t is an eigenvalue. Further, each real number is a 
nondegenerate eigenvalue of T*. The overlap of the domain of two self-adjoint exten- 
sions is the domain of T. Thus, no eigenvector of T* can be in the domains of two 
self-adjoint extensions. 

We can therefore conclude that for each simple symmetric operator T with deficiency 
indices (1, 1) whose self-adjoint extensions have purely discrete spectra, the family of 
spectra {tn{a)} of its self-adjoint extensions T{a) covers all of IR exactly once. As we 
explained above, this fact directly yields the sampling theorem of Eqs lfip?! Note that, 
since we are working in a complex Hilbert space but wish to sample real-valued and 
continuous signals, we will have to prove below that it is possible to choose the phases 
of the eigenvectors {t} such that the sampling kernel is real and continuous. 

5 Time- varying bandwidths 

5.1 Defining the time- varying bandwidth as a time- varying 
Nyquist rate 

Different choices of simple symmetric operators T lead to sampling formulas for differ- 
ent classes of signals. Let us now identify the time-varying bandwidths of these classes 
of signals, as a function of the associated family of spectra {t„(Q;)}. As discussed above, 
it is impossible to define the notion of time- varying bandwidth in terms of the signals' 
Fourier transforms. Our strategy is, therefore, to define a time-varying bandwidth 
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^^max{t) directly through a time- varying Nyquist rate, namely through a time- varying 
sampling lattice spacing function s{t). We can then define cumaxit) — l/2s{t). To this 
end, we need to show that for the spectra of the self-adjoint extensions of T the lattice 
spacing s{t) is indeed a well defined function of time. 

To see this, note that each family of lattices {tn(Q;)} is parametrizable such that the 
tnia) are differentiable, strictly monotonic functions of a, obeying, say, t„(Q; -|- 27r) = 
tn+i{a)- We will choose the C/(l)-parametrization as e'". The tn{a) are then strictly 
monotonic since, as we saw, each eigenvalue of T* is nondegenerate and no eigenvector 
can be in the domain of two self-adjoint extensions. 

Crucially, the strict monotonicity of the tn{a) implies that, for a given family 
{tn(tt)}, there is exactly one lattice interval centered around each point t. Let us 
denote the length of this interval by s{t). Thus, if 

^ ^ tn+i{a) + tnja) ^^^^ 

for some a, then we set 

s{t) := tn+i{a) - tnia) . (12) 

For each such operator T, the lattice spacing is therefore a well-defined function of 
time s{t), i.e. we arrive at a well-defined time- varying Nyquist rate. Through s{t), we 
may now define the time- varying bandwidth as 

U^maxit) (13) 

Vice versa, given Umaxit), the family {tn{a)} can be reconstructed, up to trivial 
reparametrization of a, through 

tnia) = t - J (14) 

and 

t„+i(Q;)^t+ ^ (15) 

It is, therefore, equivalent to specify the family of spectra {t„(Q;)} or to specify the 
time- varying Nyquist rate s(t), or to specify the bandwidth curve ujmaxii)- 

The conventional characterization of the bandwidth as the width of the Fourier spec- 
trum becomes applicable only when the bandwidth is locally approximately constant. 
We note that whenever the usual Fourier-based characterization of the bandwidth is 
applicable then it yields uJmax '■— l/2s and it therefore agrees with our definition. 
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5.2 Bandlimited functions form the domain of a simple sym- 
metric time operator T 

Technically, we are characterizing a class of bandlimited signals as an analytic domain 
of the "Time operator" T : T0(t) = t(f){t). In this way we were able to generalize the 
notion of bandwidth to time-varying bandwidths. There is a simple argument for why 
a class of bandlimited signals should be an analytic domain of T (i.e. a domain on 
which T can be applied an arbitrarily finitely many times). To this end, we need to 
begin with precise definitions in the case of constant bandwidth of what we mean by 
bandlimited functions and what we mean by strictly bandlimited functions: 

We define the class B^^^^ of strictly bandlimited signals with bandwidth uJmax as con- 
sisting of all those signals whose Fourier transform 0(0;) has the properties that all its 
derivatives d^(f){uj) exist and are square integrable for all n, and that there exists a 

0^0(0) < i^max SO that 

0(cu) = for all \u\ > cuo(0) (16) 
We define the scalar product of two signals through 

duM^YM^) (17) 

The closure of B^^^^^ with respect to the norm induced by EqEI yields the Hilbert 
space H^^^^ := -Bw,„„^ which we call the space of all bandlimited signals. (What we 
really mean is of course a space of physical fields that obey an ultraviolet cutoff). 

Now we observe that if a signal 0(t) is strictly bandlimited, i.e. obeys Eq lTBl then also 
the signal T(j){t) = t(j){t) is strictly bandlimited and obeys Eqll6[ This is because the 
Fourier transform of the multiplication operator T : T(j){t) = t(f){t) is the differenti- 
ation operator. Clearly the action of the differentiation operator action changes the 
support of the spectrum at most infinitesimally. Thus, if the signals Fourier transform 
has support strictly only inside the bandwidth interval, then also the derivative of the 
signals Fourier transform has support strictly only within the bandwidth interval. 

Crucially, the fact that multiplying a signal 0(t) with the time variable, 0(t) — ^ t(f){t), 
affects the bandwidth of the signal only infinitesimally, is independent of which band- 
width has been chosen. This means that, independently of the chosen bandwidth, the 
strictly bandlimited signals form an analytic domain of the time operator. 

For this reason it is clear that any proper definition of time- varying bandwidths should 
have the property that each class of signals with time-varying bandwidth forms an 
analytic domain of the time operator. 

Finally, we remark that, technically, we obtain Dt from B^^^^ by closing the opera- 
tor (self-adjoint extensions are most conveniently constructed when starting from an 
operator which is closed), and we obtain H^^^^ by closing the domain (or B^^^J). 
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5.3 The bandwidth Umits its own variabiUty 



Intuitively, consistency conditions on the possible bandwidth curves uJmax{t) are to be 
expected. This is because the bandwidth uJmax{t) poses a limit on how much signals 
can vary around time t and the bandwidth should therefore also pose a limit on how 
much it, itself, the bandwidth Umaxit), can vary around time t. Indeed, there already 
follows a simple consistency condition from its definition in EqllHI 

Starting from 

d . , 1 ds{t) 



and using 

we calculate separately 



ds{t) ds{a) da{t) 
dt da dt 



ds{a) d 



and 



da da 
da{t) ( d + t„+i(a; 



dt \da 
Thus, introducing the notation 



(19) 



{tn+i{a) - tn{a)) (20) 



(21) 



<(«) := tr{a) (22) 
da 

we obtain: 

-c^_(t) - -4a;_(t) (23) 

Eq|2ni implies, in particular, a limit which the bandwidth ujmax{t) imposes on how much 
^max(t) itself can vary: 

I du^axit)/dt I < 4 ul,,it) (24) 



5.4 The (varying) bandwidth Umits how much a bandUmited 
signal can be peaked 

We will now show that if a signal (j)(t) is strictly bandlimited by a constant Umax, and 
if (j){t) is sampled at the Nyquist rate, then the standard deviation AT of the sampled 
values is at least AT > l/Aumax- Thus, the smaller the bandwidth, the larger is the 
signals' minimum standard deviation, i.e. the less they can be peaked. A time-varying 
bandwidth means that the lower bound to the standard deviation of the signals' samples 
depends on the time t around which the signal is centered (its first moment). 
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We saw that for each Nyquist rate there is a whole family of sampling lattices. 
We will find that the standard deviation (as well as all other moments) of the sampled 
values of a strictly bandlimited signal do not depend on the lattice on which the samples 
are taken. 

The minimum standard deviation (Jmin{t) 

We have defined time-varying bandwidths as time-varying Nyquist rates. Let us now 
investigate in which sense such a time-varying bandwidth is limiting how much signals 
can vary over time. We will find that the time-varying Nyquist rate is limiting how 
much the signals can be peaked at different times. 

Consider for example a normalized signal (f){t) 

(0,0) = 1 (25) 

which is centered around the time t in the sense that t is the expectation value of the 
time operator T, i.e. its first moment: 

(0,T0) = t (26) 

A measure of how much the signal (j){t) is peaked around its expectation value is for 
example its second moment, i.e. its formal standard deviation: 

AT(0) = v'(0,T20)-(0,T0)2 (27) 

In order to explicitly evaluate EqJTTl we can insert any one of the resolutions of the 
identity 1 = ® induced by the self-adjoint extensions T(a) of T, to obtain: 

{^T{<t>)f = J2^itn{a)r tlia) - (e0(^-(«))'^-(«)) (28) 

Thus, AT(0) is the standard deviation of the discrete set of samples of the signal 0(t) 
on a sampling lattice {t„(a)}. The standard deviation of a signal's samples 0(t„(tt)) on 
the lattice {t„(a)} is clearly the same as the standard deviation of the signal's samples 
(f){tn{a')) on any other lattice {t„(a;')} of the family. 

We note, that therefore in the Hilbert space of bandlimited signals, the standard 
deviation AT{(j)) is not the usual standard deviation of a continuous function, which 
would be calculated in terms of integrals, but instead it is the standard deviation of 
the discrete set of samples of the signal, taken on a sampling lattice. 

Our claim is that for the class of signals which is characterized by some time-varying 
Nyquist rate s{t) there is a function aminit) which sets a time- varying finite lower bound 
to the standard deviation AT of the samples of signals which are centered around t: 

AT {(f)) > aminit) for all (p e Dt which obey {(j),T(f)) = t (29) 
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Concretely, this means that if a normahzed signal is sampled on a sampling lattice 
{^n(tt)} (where a is arbitrary) and if the average, or first moment, of the samples is 
t — J2n^{'tni(^)y tn{a) then we claim that the standard deviation of the samples is 
larger or equal than a finite value cr^in{t)' 



\ 



(30) 



Let us calculate the minimum standard deviation aminit) from the Nyquist rate s{t), 
i.e. from the bandwidth cUmaxit) — l/2s{t). 



Calculating (Jmin{t) 

In order to calculate (Tmin{t) for strictly bandlimitcd signals, it will be convenient to 
calculate first the minimum standard deviation curve (Jmini^) ^ larger set signals, 
namely all those in the domain DT{a) of a self-adjoint extension T{a) of T. (We recall 
that the dense set of strictly bandlimited signals is an analytic domain of T, i.e. they 
are contained in Dt)- 

In the domain DT(a) of a self-adjoint extension T{q:). there are eigenvectors t„(Q;) 
of T{a) and clearly each is centered around time t„(Q;) and has vanishing standard 
deviation. We recall that the standard deviation vanishes exactly only for eigenvectors. 
Thus, in the domain i^T(a) of a given self-adjoint extension, the minimum standard 
deviation vanishes at the eigenvalues 

^InLitnia)) = 0. (31) 

On the other hand, for all vectors in DT(a) which are centered around a point in time r 
which is in between two neighboring eigenvalues, say tn{a) < r < tn+i{cK) the standard 
deviation cannot vanish, i.e. we expect that: 

^il(r)) 7^ 0. (32) 

In order to calculate the minimum standard deviation, u^^Jni'^)-, we use that the most 
peaked signal 0^ G Dti^u) which is centered around r must be a linear combination of 
the closest two eigenvectors, namely 

0r = aitn{a) + a2tn+i{a) (33) 

The requirements {(f)r,T(f)r) ~ r and {(f)r,(l)T) — 1 determine the coefficients: oi = 
{tn+i — T)/{tn+i — tn), and 02 = (t — tn) / {tn+1 — tn)- The Standard deviation AT'(0t-) 
of the signal (pr which is maximally peaked around the time r is then easily calculated 
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and yields the minimum standard deviation cr^j„(r): 

(34) 

= ^{ct>r,T^,)-T^ (35) 

It is clear that for each choice of a the minimum standard deviation al^in is zero at 
all times {tn{a)}, since these are the eigenvalues of T(a). EqlHUl now shows that in 
the times between any two neighboring eigenvalues tn{a) and t„+i(a) the minimum 
standard deviation curve cr^/„(t) is simply a half circle with diameter equal to the 
lattice spacing s = tn+i{a) — tn(«)- 

We can now calculate the minimum standard deviation curve (Tmin{t) ^oi the for the 
vectors in Dt-, from the knowledge of the minimum standard deviation curves (T^mini^) 
of the self-adjoint extensions T{a). 

This is because Dt C -Dr(a) for all a, and therefore the minimum standard deviation 
curve (Tmin{t) for the vectors in the domain Dt is bounded from below by all the 
minimum standard deviation curves (J^mlni^) of the self-adjoint extensions. We obtain: 

(^min{t) = max af°i,(t) (36) 

Since strictly bandlimited signals are elements of Dt, this means that strictly ban- 
dlimited signals centered around time t are peaked at most as much as specified by the 
minimum standard deviation <Jmin(t)- Here, (Jmm(^) is the standard deviation of the 
discrete samples of the signal on (any arbitrary) one of the family of sampling lattices 
{tn{<y)} with the Nyquist rate s(t). 

In physical terms, this means that a finite ultraviolet cutoff of this kind, constant 
or varying, corresponds one-to-one to the existence of a constant or varying finite 
minimum to the formal uncertainty in position (or time), as mentioned above in the 
context of studies in quantum gravity and string theory, see [Tj . 

amin in the case of constant bandwidth 

In the case of a constant bandwidth, Eq|2Slis readily evaluated: 

For a constant bandwidth Umax, the lattice spacing is s = l/2umax- Thus, each curve 
^mini't) consists of half circles with radius l/4uJmax which stretch from eigenvalue to 
eigenvalue. Therefore, from Eql3(i| the minimum standard deviation curve (Tmin{t) for 
strictly bandlimited signals is a constant: 

C^minit) = (37) 

^'^raax 
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Let us confirm this result by explicitly solving the variational problem for finding the 
minimum standard deviation: 

Our aim is to minimize (AT)^ for cf) G Dt by minimizing (0, T^0) while enforcing the 
constraints (0, 0) = 1 and (0, Tcj)) = t. 

To this end, we represent T in Fourier space as T = —i d/duj on strictly bandlim- 
ited signals, i.e. on signals which in Fourier space obey, in particular, the boundary 
condition: 

4>{-^max) = = (f){uJmax) (38) 

Introducing Lagrange multipliers ki,k2, the functional to be minimized is, therefore 



dio {(9^0*)(6'^0) + h{(f)*(f) - ci) + k2{-t(j)*d^(j) - C2)} , (39) 

which yields the Euler-Lagrange equation: 

- dl4> + k,4> - ik2d^4> = (40) 

For each value of a T-expectation value t, there exists exactly one (up to phase) nor- 
malized solution which also obeys the boundary condition, namely: 

Mu^) = ^=i= cos (-^) e-- (41) 

The standard deviations AT{(j)t) of these solutions are straightforward to calculate in 
Fourier space, to obtain 

AT{(Pt) = foralH (42) 

which indeed coincides with our previous result, Eq lHTl for the minimum standard 
deviation o"mm in the class of strictly bandlimited signals of constant bandwidth Umax- 

Relation to superoscillations 

For the special case when the bandwidth is constant, the property of strictly bandlim- 
ited signals to obey a minimum standard deviation may appear to be not more and 
not less intuitive than the conventional belief that bandlimited signals are slow- varying, 
namely that bandlimited signals are varying at most as fast as their highest Fourier 
component. 

Let us emphasize, however, that the conventional belief that bandlimited signals 
cannot vary faster than their fastest Fourier component does not hold true. This has 
been pointed out in particular in ^1 E], where counterexamples have been given 
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and named "superoscillations" . In JH], Berry quotes I. Daubechi in explaining that 
phenomena related to superoscillations have also been known to occur as instabilities 
in oversampling. For related work, see also e.g. jTHl- 

On the other hand, our finding that strictly bandlimited signals obey a minimum 
standard deviation is true for all strictly bandlimited signals, without exceptions. Su- 
peroscillations are naturally accommodated in this framework in the sense that strictly 
bandlimited signals can contain arbitrarily sharp spikes, and that they can therefore 
be locally varying arbitrarily faster than their highest Fourier component, while still 
obeying the standard deviation constraint, see [T7]. 

6 Types of sampling lattices 

6.1 Unbounded and semibounded sampling lattices 

In those cases where the lattices of sampling times are ranging from — oo to +00, the 
operator T and its self-adjoint extensions are fully unbounded operators. 

Since simple symmetric operators can be semibounded, it is also possible to choose 
sampling lattices which are merely semibounded, i.e. which range from —00 to some 
finite largest sampling time t^, or of course from some finite sampling time to +00. 
However, we found above, on functional analytical grounds, that the family of spectra 
always covers all of IR, i.e. the sampling formula always allows to reconstruct for 
all t e IR. It appears to be paradoxical that signals can be sampled on a semibounded 
lattice but recovered for all t G IR. Nevertheless, it is a fact that simple symmetric 
operators can be semi-bounded and that in such cases can still be reconstructed 
for alH G IR from its values sampled on a semibounded lattice. 

The resolution to the paradox is that the bandwidth curve Wmaxit) rapidly decays 
after the last sampling point t^. This means that the class of signals which are sampled 
in such cases is a class of signals which quickly become constant in the interval from 

to infinity. 

To see this, assume that the spectrum of the self-adjoint extension T(0) is, say, 
bounded from above by possessing a largest eigenvalue tN- We recall that each real 
number appears as an eigenvalue in one of the self-adjoint extensions of T. This means 
that as a runs through the family of self-adjoint extensions, the highest eigenvalue 
t]\f{a) of their respective spectra moves from tAr(O) = t^r to infinity, t^^a) +00 for 
a 27r, while the next to highest eigenvalue behaves as tM^i{a) — > t]^. 

Thus, while a runs through the family of self-adjoint extensions, the center t = 
(tAr(Q;) — t7v-i(tt))/2 of the interval between the last two eigenvalues moves to infinity, 
and thereby the spacing s(t) = tAr(Q;) — tAr_i(a;) of this interval diverges. Thus, the 
bandwidth Umaxit) = l/2s(t) indeed decays like uJmaxit) — > 1/t between the last point 
of the sampling lattice and infinity. 

We will show in Sec lfi.HI that this phenomena can be used to sample signals on 
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lattices which are not actually finite lattices, but which are effectively finite lattices for 
all practical purposes. 



6.2 No finite sampling lattices 

Before we discuss sampling on effectively finite sampling lattices, let us convince our- 
selves that there are no actually finite lattices of sampling times: 

To this end, we recall that we are defining the class of signals to be sampled as 
the representation of the multiplication operator T : T(f){t) = t(j){t) on the space of 
functions for which the scalar product 

(01,02) :=E'/'I(^n(«)) UUa)) (43) 

n 

is independent of a, where {t„(a)} are the spectra of a family of self-adjoint extensions. 
For example, the Shannon sampling formula is now the expansion of signals 

m = Y.<^^:^t\t) (44) 

n 

with coefficients 

ci") = 0(t„(a)) (45) 

in the Hilbert basis of functions: 

27r [t - tn{a)) Umax 

The basis functions bl^\t) are orthogonal with respect to the scalar product given by 
EqliHl Here, {t„(a)} is the family of lattices with the equidistant lattice points defined 
as in Eq|21 through t„(a) = a + n/2umax- We recall that for vectors in the Hilbert 
space, the scalar product EqES] yields the same result for any arbitrary choice of a. 

Here, and also in the case of generic nonequidistant sampling lattices {t„(a)}, the 
a-independence of the scalar product implies the existence of isometries between the 
/^-Hilbert spaces over each of the lattices of the family {tn{a)}- The matrix elements 
of these isometries then form the sampling kernel G. 

This means that the signals are solutions to the equation 

<P{tm{a')) = G{tm{a')Mc^)) <P{tn{c^)), (47) 

n 

for all a and a' . 

In addition, the signals (j) form a representation of T, i.e. they also obey 

tm{a') (t){tm{a')) = J2 G(t„(a'),t„(a)) t„(a) 0(t„(a)) (48) 
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for all a and a'. Assume now that the family {^^(q;)} consists of finite lattices, say 
with lattice sites. The sampling kernel G is then a x A^- dimensional matrix with 
matrix elements «) := G{trn{(^')-,tn{oi))- In this case, Eqs J47|481 would yield 

the following equation for the sampling kernel matrix G = Gnm{o-' ■, o)'- 



It is clear that this equation cannot be solved in finite dimensions, since the adjoint 
action of any isometric matrix G would be spectrum preserving. This simply refiects the 
functional analytic fact that finite dimensional symmetric operators are automatically 
self-adjoint, i.e. that only unbounded operators can be simple symmetric. 

6.3 Effectively finite sampling lattices 

In practical applications in engineering it is only possible to use finite lattices to sample 
actual signals, i.e. the theoretically infinite lattices are cut off. 

The generalized sampling theorem allows us to perform the cutoff of lattices in a sim- 
ple and controlled way. Namely, outside the finite time span, say [tzo5^«p] in which 
actual sampling takes place, we are free to choose arbitrarily large lattice spacings. By 
choosing very large lattice spacings outside the interval \tio,tup\ we achieve that the 
bandwidth quickly decays to zero outside the interval. More precisely, the bandwidth 
decays like l/\t — tup\ and l/\t — tio\ to the left and right of the boundaries of the 
interval [t/o,^«p], meaning that we are sampling signals which quickly become constant 
outside the interval. Since the signals are in the Hilbert space and are therefore square 
summable, they indeed decay to the constant zero. This shows that while the un- 
boundedness of the operator T is a mathematical necessity this does not affect in any 
essential way the sampling properties of the signals (or fields) that are based on finite 
sets of samples. After we develop the details of the sampling theorem we will come 
back to effectively finite sampling lattices in Sec|H] 

7 Calculating the sampling kernel 

7.1 Embedding a self-adjoint operator into a U{1)— family of 
self-adjoint operators 

We begin by specifying a lattice of sampling times on which signals are to be 
sampled and we define the self-adjoint operator T(0) which has this lattice as its 
spectrum. 




(49) 
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Our strategy is to define a self-adjoint operator T(0) with the lattice {tn} as its 
spectrum, then to use the data {i^} to specify a simple symmetric restriction T of 
T(0), and finally to construct its U{1)- family of self-adjoint extensions T{a). The 
scalar product of their eigenvectors yields the desired sampling kernel G{t,tn)- 

To this end, consider the self-adjoint operator 

T{0):^J2intn(S)K (50) 

n 

on its domain Dr(o)- Here, t„ denotes the eigenvector to the eigenvalue tn, and t* 
denotes the dual vector. The Cayley transform of T{0) is a unitary operator 

m-Er^^n^K (51) 

which is defined on the entire Hilbert space H = Dt(o). We restrict U{0) to a simple 
isometric operator 

S = UmHec.+ (52) 

by taking from its domain one dimension spanned by a normalized vector v_|_ G H. 
The inverse Cayley transform then yields a simple symmetric operator T on Dt 

Dt = {0 e Dt(o} I (v+, (r(0) + i)cj>) = o} (53) 

iff the so-defined Dt is dense and contains no eigenvectors. In the {t„}-basis, the 
condition (v+, (T'(O) -|- i)4>) = reads Z^n(v+, t„) {tn + i) (t„, 0) = 0. Thus, for Dt to 
be dense, v+ must be chosen such that 

E I (v+, tn) {tn + 01' = divergent (54) 

n 

since otherwise all of Dt would be orthogonal to w := Z^n(tn + i) (v+,tn) tn, which 
would then be a vector in H. 

Further, the condition that there are no eigenvectors in Dt implies that v_|_ must 
be chosen such that (v+,t„) ^ for all n. The second deficiency space is of course 
spanned by v_ = ^7(0) v_|_. Now a t/(l)-family of self-adjoint extensions T{a) of T is 
obtained as the inverse Cayley transforms of the unitary extensions U {a) of 5": 

U{a) = (l - v_ vl + e'" v_ (g) vl) U{0) (55) 

The choice of a spectrum {tn} = {^n(0)} and a deficiency vector v+ thus determines 
a family of self-adjoints T{a) and unitaries U{a) with their eigenbases {t„(Q;)} (up 
to phases). The functions tn{(x) are strictly monotonic since each cigcnvahic of T* is 
nondegenerate and no eigenvector can be in the domain of two self-adjoint extensions. 
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7.2 Calculating the ^-dependence of the sampling kernel 

The aim is to calculate the sampling kernel, i.e. the scalar products of the eigenvectors 
of the unitary and self-adjoint extensions. To this end, let tr{a) be a normalized 
eigenvector of the operators U{a) and T{a). Then: 

544^ (tn,t.(a)) = (t„,C/(a)t,(a)) 

= (t„, (1 + (e*" - 1) v_ ® U{0) tr{a 



tn + i 



(t„, tr{a)) + (e'" - 1) (t„, v_) (v+, tr{a)) 



Thus, 



where we defined 



and 



N{tr{a)) := ^-r^^ (e^- - 1) (v+, t,(a)) (57) 



/„ := (tn + t) (t„,v_) (58) 



In this way we have separated the t,,(a)-dependcncc of the sampling kernel in the form 
of N{tr{a)) from its n-dependence through the coefficients /„. It is now straightforward 
to calculate N{tr{a)), i.e. the tr(«)-dependence of the sampling kernel by using the 
normalization condition: 

1 = Z^(tr(Qi),t„j (t„,t^(a)) = 2^ _ (59) 

n n VrKO^) '"nj 

Thus, for all t in the family {^^(0;)}, i-e. for all t e IR, there holds 

7.3 Calculating the n-dependence of the sampling kernel 

In order to calculate the coefficients /„, we begin by noting that 

(t„,v_) = (W,t/(0)v+) 



= J2 (t„,^7(0)tj(t^,v+) 

m 



19 



so that, in addition to EqEHl the /„ also obey 

fn = {tn - i) (t„, v+) (61) 

We will now show that oc /*/„, i.e., more precisely, that dtn{(y)n/ da\a=o = f^fnl'^- 
To this end, inserting a resolution of the identity in terms of the {t„}-eigenbasis 
into EqlHZI yields 

Ar(t,(«)) = M^(e--i)^ (v+,t„) (t„,t,(a)) (62) 
Inserting EqlSHland Eq ]611 we obtain 

N{tM)) = M^(e- - 1) E , ^^^^^^"^^ ^ (63) 



Thus: 

2i ^ fUn {tr{a)+i) 
e^"-l ^ {tn + i){tr{a)-tn) 



(64) 



Let us calculate separately the imaginary and the real part of Eq]64l The imaginary 
part yields: 

1 , 1 - ^ f*fn (^(")+^)(^"~^)~(^(")~^)(^" + ^) (65) 



e- _ 1 e— - 1 2z V itl + - tn) 

and thus, 

Using Eq l611 we see that Eq 1661 simply expresses the normalization condition for v+. 
On the other hand, the real part of EqJ64l yields: 

e--i e— -1 2^ 4^^"^"M^(t„ + z)(t(a)-tO^(t„-0(t(a)-t„)y'' ^ ^ 

1 p . (^(«) + - ^) + (^(«) - ^)(^« + ^) ..o^ 
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and thus, using 



gia/2 _j_ g-ia/2 

cot(a/2) = i —nr 



i 



2 - 



we obtain: 



Since — ^ for a ^ 0, EqlHni finally yields: 

dtr{a) f*fr 



da 



1 0=0 



(70) 



7.4 The deficiency vector v+ is parametrized by the coeffi- 
cients fn 

Through EqlHH the deficiency vector v+ is parametrized by the coefficients /„. Thus, 
the constraints on the choice of v_|_ are in one-to-one correspondence to constraints on 
the choice of coefficients /„, namely 

f*.fn 

fn^O forall n, ^ " " = 1 and Y.fnfn = divergent . (71) 

Clearly, such /„ only exist if |t„| ^ 00 for n — > 00 and/or for n ^ —00. This refiects 
the functional analytical fact that only semibounded or fully unbounded symmetric 
operators can have nonzero deficiency indices. 

We can now use EqEOlto express the conditions on v+, as conditions on the coef- 
ficients {t^}. We obtain as the constraints on the data and {t^}'- 

tn+i > tn , and > for all n, (72) 

— > 00 for n ^ +00 and/or for n ^ —00 (73) 

and 

Et1^= finite, divergent. (74) 

7.5 Realness and differentiability of the sampling kernel 

While we are working in a complex Hilbert space, the aim in engineering would be to 
sample real-valued signals. In the case of physical fields, we may want to work with 
charged, i.e., complex- valued scalar fields, but we will certainly also want to be able to 
work with real-valued scalar fields. 
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We therefore have to deal with the fact that the representation of vectors of the 
Hilbert space as functions = (t, (p) a priori depends on the arbitrary choice of the 
phases of the vectors t. Correspondingly, also the sampling kernel (t,t„) depends on 
the choice of these phases. For the sampling of real-valued signals we therefore need 
to show that it is possible to choose the phases of the vectors t such that if a signal 
is real-valued on one sampling lattice {tn}, then it is real-valued for all t. This means 
that we need to show that it is possible to choose the phases such that the sampling 
kernel is real-valued. 

In order to derive the appropriate choice of phases, we reconsider the sampling 
kernel Eg 1561 

(tn,t) = 

which holds for all t tn, while, of course. 

We observe that the denominator in EqlTSlgoes linearly to zero whenever t approaches 
one of the tn of the lattice. Nevertheless, the absolute value of the sampling kernel does 
not diverge since, as follows from Eg 1601 simultaneously also the numerator linearly 
approaches zero as t approaches tn- 

In particular, whenever t passes through a t„ the denominator in EqESl changes 
sign. Thus, in order for the sampling kernel to be continuous, it is necessary that also 
N(t) changes sign at the points tn, i.e. as it passes through a zero. 

Indeed, the definition of N{t) in Eg IFTI shows that it is possible to choose the phases 
of the vectors t 7^ t„ such that 

N{tr{a)) = {-iy\N{tr{a))\ . (77) 

This choice of phases ensures that the sampling kernel is real and continuous for all 
t 7^ tn, and that, for all n, both limits of the sampling kernel (t„,t) for t — > t„, from 
above and below, coincide. However, in order to obtain that the sampling kernel is 
continuous at the times t„, we still need to ensure that those limits also coincide with 
the value (t„,t„) = 1 set by the normalization condition Eg]76l Therefore, since N{t) 
changes sign at each tn, also the /„ must have alternating signs: 

/„ = (-l)'^ \fn\ (78) 

We still have the freedom to choose the phases of the vectors t„. Indeed, using EgEH 
we can choose the phases of these vectors so that Eg I7H1 holds. 

We remark that the choice of phases to obtain the sampling kernel is unigue up to a 
global phase: The sampling kernel is a scalar product of the vectors t and therefore the 
same real and differentiable sampling kernel is obtained if all vectors t are multiplied 
by some constant overall phase c. Indeed, the definitions of N(t) and show that 



fn N{t) 
t tn. 



1 . 



(75) 



(76) 
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this amounts to multiplying N{t) with c and multiplying /„ with c*, so that the phase 
cancels out in the expression for the sampling kernel. 

We have therefore shown that and how the phases of the vectors t can be chosen so 
that the sampling kernel is real and continuous, to obtain: 

(t„,t^(a)) = (-1) . X _ (79) 

The sign factor can be written, equivalently, as 

where z{tr{a),tn) is the number of sampling points tj in the interval [tr{a),tn\ if tr{a) < 
tn or the number of ti in the interval t] if t„ < tr{a), and z = if t„ = tr{(y). 

Since each t G IR occurs in the family of lattices {tr(a)} we can now express the 
sampling kernel in the general form: 

(t„,t) = (-i)-(*A) \fn\\m\ ^g^^ 

Finally, using or results Eqs l70l57| which express and |A^(t)| in terms of the data 
{tn} and {t^} we obtain: 

t tn \ m i.^ tm) J 

7.6 Calculating the full sampling kernel G{t,t'), which yields 
all sampling lattices and Umax{t) for all t 

So-far we have only calculated the sampling kernel (t,t„) for one sampling lattice 
{tn} = {tn{0)} of the family of sampling lattices. It is straightforward to calculate the 
general sampling kernel: 

(t,,(a),t,(a')) = E (tr(a),t„) (t„,t,(a')) 

n 

(-1Y+' tL 



or, in simplified notation, for arbitrary times t, r G IR: 



(t tn)(r tn)^' (Em {t-L?) (r-L)O 
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The significance of this equation is that, given one samphng lattice {t„} = {tn(0)} and 
its derivative {t^} = {^^n(0)}, Eg IH^ allows us reconstruct the entire family {tn{a)} 
of sampling lattices, up to reparametrization of a. This is because each of the sampling 
lattices consists of the eigenvalues of orthogonal vectors, i.e. of vectors whose scalar 
product vanishes: 

For any to ^ IR there exists a self-adjoint extension T{a) and some index n so that 
to = The entire sampling lattice {tn{a)} then consists of to and all t for which 

(to,t) =0: 

{tn(«)} = {to} U {t I (to,t) = 0} (84) 

To this end, for generic data {t„} and {t'^} the zeros of the general sampling kernel in 
EqlHSl can at least be found numerically. Thus, all lattices of the family {t„(a)} can 
be calculated, in principle, from the data {tn}, {t^} which enter the sampling theorem. 

In turn, once all sampling lattices are known, also the curve LJmaxit) can be con- 
structed explicitly for all t, as explained above. 



7.7 Other sampling lattices 

In the new sampling theorem, the class of signals to be sampled is chosen by specifying 
one samphng lattice {tn} and its derivative, {t'^} and by thereby specifying the time- 
varying bandwidth curve ujmax{t)- The sampling theorem allows us to recover the 
signals of this class from their values at any one of the sampling lattices from the 
family of sampling lattices which is determined by the data {t„} and {t^}. 

In fact, the same class of signals, with bandwidth specified by the data {t„} and 
{t'n}, can also be sampled on a large variety of other sampling lattices, by using a 
suitably modified sampling formula. 

To this end, assume that {t„} is any sampling lattice for which there exists a linear 
transformation K, such that: 



Applied to 
this yields 



Ksr (tr)tm) — Sgm (§5) 

r 

(t,,t^) {tm,<P) = (tr,0) (86) 

m 

(t„,0) = ^n,r {tr,tm) (tm, 0) 

m,r 

= ^n^r (tr,0) 



Inserting this result into (t,0) = J2n (tjtn) (t„,0), we obtain the sampling formula 
for sampling on the lattice {t„}: 

{t,(P) = Y (t,tn) Kn,r (t., 0) (87) 
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8 Summary of the new sampling method 



The new sampling method allows us to continuously increase and lower the sampling 
rate over time in order to sample signals with a correspondingly time-varying band- 
width. For engineering purposes, this allows one to economize on data collection when- 
ever the bandwidth is temporarily low and to increase the sampling rate whenever the 
bandwidth is temporarily high. This is the basis of the data compression method 
patented in 0. 

The mathematical background is based on the observation that the action of the 
"time" operator T : T(f){t) = t(f){t) can only infinitesimally change the spectrum of 
bandlimited signals (because it acts as differentiation in Fourier space), independently 
of the chosen bandwidth. We were able, therefore, to identify classes of signals with 
time-varying bandwidth as analytic domains of simple symmetric multiplication oper- 
ators T : — >• tcpit). (In the case of fields in spacetime, we may think, equivalently, 
of T position coordinate.) 

In particular, we found that, as in the case of constant bandwidth, also in the case 
of a time- varying bandwidth, there exists a family of lattices {t„(a)} of sampling times 
at which the signals can be measured and recovered for all t. We found that the lattice 
spacings of these sampling lattices, together, form a well-defined function of time s{t). 
This yielded a precise definition of time-varying bandwidths as time-varying Nyquist 
rates, i.e. as a time-dependence of the spacings s{t) of sampling lattices. We defined 

UJmax{t) = l/2s(t). 

We also discussed that, contrary to intuition, a finite bandwidth does not imply that 
signals vary "at most as fast as their highest Fourier component". We found that 
instead a true property of finite bandwidth signals is that they obey a lower bound to 
how much they can be peaked, in the sense that the standard deviation of the samples 
(j){tn{c()) (for arbitrary fixed a) of strictly bandlimited signals is bounded from below. 

More generally, we found that classes of signals with time- varying bandwidths obey 
a time-varying lower bound on their standard deviation, i.e. they consist of signals 
which, if centered around time t, have standard deviations bounded from below by a 
positive time-dependent function (Jmaxit)- 

8.1 How to use the sampling theorem 

Concretely, the sampling theorem requires as input data first the choice of one sampling 
lattice {tn}- The choice of one sampling lattice largely determines the Nyquist rate and 
therefore the time- varying bandwidth uJmax{t)- This is because by giving one sampling 
lattice the Nyquist rate s{t), i.e. the lattice spacing, is defined already for a discrete 
set of intervals. This means that the bandwidth is already defined at a discrete set of 
times, namely at the times r„ = (t„ + t„+i)/2, i.e. at the center of each lattice interval. 



25 



through the inverse of the spacing: 

We also found that the bandwidth cannot vary much between these points: 

d 



7, ^max(t^ 

dt 



< 4 UJmaxity. (89) 



The samphng theorem requires a further set of input data, {t^}, to fine-tune the 
bandwidth curve by specifying the derivative of the bandwidth at those points, namely 
through: 

— I ) (t \ — —At) r-r "l^ ^^"+1 ~ ^n) (Qrs\ 
J, ^max\'n) ^ ^maxy'n) /,/ \ W'-'/ 

at ytn+l + '^n) 

We showed that the bandwidth curve uimax{t) thereby becomes determined for all t. 
We recall that by embedding the sampling lattice {tn\ as the lattice {tn(0)} into 
the family of sampling lattices {tn(a)} which is determined by (and determines) the 
bandwidth Umaxit), the auxihary data can be identified as the derivative of the lattice: 
{4} = {d/dt t„(a)U=o}. 

After choosing sampling times {t„}, which already largely determine the bandwidth, 
and auxiliary data {t^}, which fine-tune and fix the bandwidth curve oJmaxif) for all t, 
the sampling theorem then allows to recover all signals of the class with this bandwidth 
through the sampling formula: 

(t>{t) = E G{t,tn) (p{tn) (91) 
n 

with the sampling kernel G{t,tn): 

G{t,Q = (_i)^(M.)^ 77^] (92) 



The number of eigenvalues t„ in a given interval defines what we call the local den- 
sity of degrees of freedom and variable bandwidth. Taking samples at the so-defined 
Nyquist rate expands functions in terms of an orthonormal basis of functions. The 
reconstruction is, therefore, stable in the sense that small perturbations in the sample 
values lead to small perturbations in the reconstructed function. Stability of the re- 
construction is important since sample values are generally subject to noise and/or to 
quantum uncertainties. For the case of conventional Shannon sampling it is known that 
reconstruction from samples taken at a lower rate than the Nyquist rate is possible, in 
principle, but must be unstable, as Landau proved in [TH]. 

Here too, we may choose to use the data {tn},{in} merely in order to fix the 
bandwidth curve oJmaxif) and then decide to sample the so-specified class of signals on 
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a different sampling lattice. Indeed, we found that the same class of signals can be 
sampled also on any other lattice, say {in}, for which Kr^s = {G{tr,ts))r,s is invertible. 
The sampling formula then reads: 

= Y: Git, t„) K-), <l>iim) (93) 

n,m 

The spacing of these sampling lattices {t^} may be larger than the spacing of the 
eigenvalues of the self-adjoint extensions of T. Since we are covering conventional 
Shannon sampling as a special case, see Sec l8.2[ it is to be expected that reconstruction 
is necessarily unstable if the samples are taken at a lower rate than the density of the 
degrees of freedom. It should be interesting to prove this in generality. 

Finally, let us briefly address the case of engineering applications, where all sampling 
lattices are finite. 

We recall from Sec J6.2l that finite sampling lattices do not occur in the theoretical 
framework. Technically, the reason is that the underlying time operator T can only be 
simple symmetric if it is unbounded. Thus, in principle, in Eqs l91l9^ the indices n and 
m which label the sampling times t„ must run through a countably infinite index set, 
e.g. the natural numbers or the integers. Therefore, a truncation error appears to be 
unavoidable. 

Our generalized sampling theorem allows us, however, to freely choose the sampling 
points {tn}, as long as the sequence {tn} is divergent. This freedom can be used to 
sample on effectively finite sampling lattices, say {tn}nZi- 

To this end, we extend the finite sampling lattice {tn}nZi on which we wish to 
sample to an infinite sampling lattice by adding an infinite number of further sampling 
points. These additional sampling points we are free to choose for example at very 
very large t, where the signals have negligible support. Thus, these sampling points 
do not contribute to either of the sums in Eqs ]91l9^ In practise, we can therefore 
leave Eqs l91l9^ unchanged and we simply restrict the sums over n and m to the finite 
sampling lattice {tn}nZi ■ 

As discussed in Sec J6.31 the class of signals which is being sampled in this case has 
a bandwidth curve Umaxit) which decays as 1/t outside the interval [ti,tAr] in which 
the samples are taken. 

8.2 Recovering conventional Shannon sampling as a special 
case 

The standard Shannon Sampling Formula is to arise as a special case. The class of 

signals with constant bandwidth Umax is characterized by the family of lattices 

with tn(«) = a + n/2ujmax- Choose the sampling lattice = {tniS^)}- Its derivative 
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{t'^} is t'^ = 1. Thus, EqinSl yields the samphng formula: 

\m=-oo v'' ''mj / n=-oo ^ ^ri 

Since this equation must coincide with the Shannon formula, Eq|21 we can conclude 
that 

7r2 g 1 



sin2(7rx) (x - n)2 

which provides an example of a series expansion deduced with the new method. Note 



that setting x = 1/2 we obtain the expansion n = JYl,n=-oQiP' ~ 1/2) 



9 On the relation to coherent states and wavelets 

9.1 One - and two parameter over complete, continuous and 
normalizable resolutions of the identity 

We are here obtaining resolutions of the identity in terms of an overcomplete and 
continuously parametrized set of normalizable vectors, 

1 r'^'^ 

ZTT JO „ 

1 f+°° , da , , , 

In this sense, our expansion is similar to coherent states and wavelet expansions. The 
main difference is that while we arrive at a one-parameter resolution of the identity in 
the Hilbert space, coherent states and wavelets, see e.g. [in], are typically characterized, 
analogously, by two-parameter resolutions of the identity. In this context, it will be 
interesting to follow the analogous construction to the one here, but for time-operators 
with higher deficiency indices. 

The resolution of the identity obtained in Eq|M] implies that we are in fact dealing 
with a reproducing kernel Hilbert space. Namely, the signals obey 

(t,0) = - dr^ir) (t,r) (r, 0) (95) 

so that, with the reproducing kernel 

i?(t,r) := (t,r) ^(r) (96) 

we obtain: 

1 r+00 

ct>{t) = — dr Rit,T) <f>{T) (97) 

ZTT J-00 
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9.2 Calculation of the reproducing kernel 

In Sec J7.fjl we calculated the sampling kernel (t, r) for all t, t in terms of the input data 
{tn}, {t'n}- In order to determine the reproducing kernel it remains to calculate da/dt. 
To this end, we solve EqEHlfor a, by using 



so that 



where 



We can now calculate 



from 



and 



Thus, we obtain 



cot" 



1 , iz + 1 
- m . 

2 iz-1' 



a 



i In 



da 
It 

da 
dy 



^yjt) + 1 

zy{t) - 1 
ttn + l 

da dy 
dy dt 



dy 
dt 



l+y^ 
1 



(t - tnf 



da 
'dt 



2 V f 



1 



2 ' 



(t-t„)(t2+l), 

which finally yields a fully explicit expression for the reproducing kernel: 



R{t, 



T 



2 t' 



(t-tn){T-U) 



(t-ts)2 



(98) 

(99) 

(100) 
(101) 
(102) 
(103) 

(104) 



(105) 



10 Outlook 

In the literature, studies of quantum gravity and string theory, see [T], have suggested 
a natural ultraviolet cutoff in the form of a finite lower bound to the resolution of 
distances with a corresponding modification of the uncertainty relations. As discussed 
in Sec ]5.4l theories that are ultraviolet cutoff in the sampling theoretic sense concretely 
realize such generalized uncertainty relations. A sampling theoretic ultraviolet cutoff at 
the Planck scale also corresponds to a finite density of degrees of freedom for physical 
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fields in the very same sense that bandhmited music signals possess a finite density of 
degrees of freedom and a correspondingly finite information density. 

This leads to the prospect that sampling theory could allow us to formulate physical 
theories which combine the virtues of lattice theories with the virtues of continuum 
theories. Fields, such as those of quantum field theoretic path integrals, could be 
defined over a continuous spacetime manifold or, entirely equivalently, over any lattice 
which is sufficiently densely spaced. A theory's ultraviolet finiteness would be displayed 
when formulated on one of the lattices, while the same theory's preservation of external 
symmetries would be displayed when written as a continuum theory. 

Here, we proved sampling results for the one-dimensional preparation for 

the generalization of Shannon sampling theory to generic curved spacetimes which has 
been outlined in [Hj and jS] • To this end, we generalized sampling theory using the func- 
tional analysis of simple symmetric operators which have deficiency indices (1,1) and 
whose self-adjoint extensions have purely discrete spectra. It should be straightforward 
to apply these functional analytic methods to simple symmetric operators with higher 
deficiency indices, thereby leading to sampling methods for vector- and spinor-valued 
fields. The very same methods should be useful for developing the details of the general 
sampling theory of fields on n-dimensional curved spacetimes. 

It is encouraging that the mathematical language of the sampling theoretical Planck 
scale cutoff, namely functional analysis, is also the language of quantum theory. 

Acknowledgement: The author is happy to thank John Klauder for useful criticisms. 
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